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Abstract 

Maximum Likelihood (ML) algorithms, for the joint estimation of synchronization impairments and channel in 
Multiple Input Multiple Output-Orthogonal Frequency Division Multiplexing (MIMO-OFDM) system, are investigated 
in this work. A system model that takes into account the effects of carrier frequency offset, sampling frequency offset, 
symbol timing error, and channel impulse response is formulated. Cramer-Rao Lower Bounds for the estimation of 
continuous parameters are derived, which show the coupling effect among different impairments and the significance 
of the joint estimation. We propose an ML algorithm for the estimation of synchronization impairments and channel 
together, using grid search method. To reduce the complexity of the joint grid search in ML algorithm, a Modified 
ML (MML) algorithm with multiple one-dimensional searches is also proposed. Further, a Stage-wise ML (SML) 
algorithm using existing algorithms, which estimate fewer number of parameters, is also proposed. Performance of 
the estimation algorithms is studied through numerical simulations and it is found that the proposed ML and MML 
algorithms exhibit better performance than SML algorithm. 

Index Terms 

MEMO, OFDM, Synchronization, Channel Impulse Response, Carrier Frequency Offset, Sampling Frequency 
Offset, Symbol Timing Error, Cramer-Rao Lower Bound. 

I. Introduction 

The integration of Multiple Input Multiple Output (MIMO) and Orthogonal Frequency Division Multiplexing 
(OFDM) techniques has become a preferred solution for the high rate wireless technologies due to its high spectral 
efficiency, robustness to frequency selective fading, increased diversity gain, and enhanced system capacity. The 
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main drawback of OFDM based systems is their susceptibility to synchronization impairments such as Carrier 
Frequency Offset (CFO), Sampling Frequency Offset (SFO) and Symbol Timing Error (STE) fl], EJ. 

So far, most studies on OFDM systems have considered synchronization impairments and channel separately 0- 
IS). Maximum Likelihood (ML) estimation algorithms for joint estimation of CFO, channel, and STE of an 
Orthogonal Frequency Division Multiplexing Access (OFDMA) system and Single Input Single Output (SISO)- 
OFDM system are proposed in (31 and ||5|, respectively, but SFO is assumed to be zero. Similarly, a joint ML 
time frequency synchronization and channel estimation algorithm for MIMO-OFDM systems has been proposed 
in H and Q, without considering the effect of SFO. In ID, a pilot-aided joint Channel Impulse Response (CIR), 
CFO, and SFO estimation scheme has been proposed for MIMO-OFDM systems, assuming a perfect start of frame 
detection. An ML estimator of SFO and CFO for an OFDM system is developed in ||9l, but STE is assumed to be 
zero with perfect channel knowledge. Joint ML estimators for SFO and channel for OFDM systems are proposed 
in Is], without considering the effect of CFO and STE. 

In this paper, we propose ML algorithms for the joint estimation of SFO, CFO, STE, and channel in MIMO- 
OFDM system. An ML algorithm is proposed, where the multi-dimensional optimization problem for estimating the 
parameters is reduced to a two-dimensional and a one-dimensional grid search. To reduce the complexity further, a 
Modified ML (MML) algorithm, which involves only multiple one-dimensional searches, is also proposed. Further, 
a Stage-wise ML (SML) algorithm using existing algorithms, which estimate fewer number of parameters, is also 
proposed. Cramer-Rao Lower Bounds (CRLB) for the estimation of continuous parameters are derived, which show 
the coupling effect among different impairments and the significance of the joint estimation algorithms. Some results 
of this paper are also presented in [ 101 . 

II. SYSTEM MODEL 

Consider a MIMO-OFDM system with Nt transmit antennas and iVfl receive antennas using Quaternary Phase 
Shift Keying (QPSK) Modulation and N subcarriers per antenna. The input bit stream is first multiplexed in 
space and time before being grouped by the serial-to-parallel converter. After Inverse Fast Fourier Transform 
(IFFT) operation, cyclic prefix (CP) insertion and digital-to-analog conversion, the transmitted signal from m*^ 
radio frequency (RF) transmit antenna undergoes fading by the channel before reaching the v^^ RF receive antenna. 
Let T be the sampling time of the sampling frequency oscillator at the transmitter Then the subcarrier spacing is 
given by \/{NT). The CIR between u"^ transmit and v^^ receive antenna is /i„,t,(r) = 'Yld=o ^ hu.v,i5{T — ti). 

Notations: Upper case bold letters denote matrices and lower case bold letters denote column vectors. A denotes the estimate of A. 5(i(A) 
and ^(A) denote the real and imaginary parts of the elements of A, respectively. 1 and represent the all-one and all-zero column vector, 
respectively. represents the partial derivative of A with respect to e. Ijv/ denotes an M X Af identity matrix. A*, A-^, A^, and A^ 

denote complex conjugate, transpose, conjugate transpose, and pseudo-inverse of A, respectively. [A]m,n denotes the (m, n)"^ element of A. 
[A]xxP represents A with K rows and P columns, (g) and o represent Kronecker product and Hadamard product, respectively and || x ||p 
denotes /p-norm of x. diag[ii\ represents a diagonal matrix having the elements of x as diagonal elements and diag[A.] denotes a column 
vector with diagonal elements of A as its elements. Tr(A) represents sum of the diagonal elements of A and ppm represents parts per million. 
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where hu,v,i denotes the channel coefficient, t/ denotes the l^^ channel path delay (t; — IT), and Lu^v denotes the 
length of CIR, for u = 1, 2, 3, ... , Nt and u 1, 2, 3, . . . , Nr. 

At each receive antenna, a superposition of faded signals from all transmit antennas together with noise is 
received. Frequency differences between RF oscillators used in the MIMO-OFDM transmitter and receiver, and 
channel induced Doppler shifts cause a net CFO of A/c in the received signal, where fc is the operating radio 
carrier frequency of the RF oscillator. Furthermore, at the receiver, the received signal is sampled at T where 
T' and AT = T' - T, which results in SFO. The normalized CFO is e = AfcNT and the normalized SFO is 
f] = AT/T. The frame arrival detection in OFDM based systems is done using different correlation methods ifTTl - 
lfT3l . The main drawback of these methods is that the correlation functions do not produce a sharp peak at the 
arrival of the frame, which results in difficulty to find the fine frame time arrival instant, causing STE ||6l. Let the 
STE, after start of frame detection, be given by integer number of samples 9T, where 9 is the normalized STE. 
The fractional part of the timing error is incorporated into the CIR lfT4l . ifTSll . We assume all receive antennas 
experience common synchronization impairments in a single user MIMO-OFDM system |l6l, Q. 

We consider the estimation of impairments during the training block. As the training blocks are usually preceded 
by long CP in practical applications, we assume that the length of CP is greater than {L^ + Omax)^ where = 
max{L„ „} and 6max is the maximum STE The received signal at the u*'' RF receive antenna is given by 

r^(n) = exp(j27rA/cnT ) ^ ^ hu,v,iSu{nT - OT ~ ti) + Wy{n). 

ti=l 1=0 

where Su{n,) is the signal transmitted by the u^^ transmit antenna and Wv{n) is the complex additive Gaussian 
noise at the v^^ receive antenna with mean zero and variance cr^. We have, 

SuinT -OT- n) = s„(n(r + AT) ^ OT - n) = -0- l)T) 

exp(j27rA/cnT') = exp(j27r£n(T + AT) /NT) = exp(j27re^n/iV). 

where — e{l + rj), and = n(l + 77). Thus, 

r^{n) = exp{j2TTe,jn/N)^ ^ hu,v,iSu{nrj - 9 - I) + Wy{n). (1) 
ti=i 1=0 

From ([TJ, it can be observed that there is coupling between the parameters e, 77, 6 and hu^v The Channel 
Frequency Response (CFR) is expressed as hu.v{k) = J2f=a ^ hu^v,i exp{~j2Trkl/N). Let Xu{n) ~ Sti(n,,). The 
frequency domain signal samples are i„(fc) — '^^^^ Xu{n) exp{—j2'Kkn/N). After removal of CP, r„(n) in ([T]i 
can be expressed as, 

rvin) = ^Mj^T^pn/N) ^ ^ expij27:n^k/N) exp{-j2Trek/N)K^,ik)i,,{k) + (2) 

u=i k=a 

where the initial offsets due to CP are assumed to be zero. Let matrices [FiJjvxAt and [F2]ivxL„ be defined as 

rp exp{j2TTk{n{l + 7j))/N) 

[tl[V)\nM = , (3) 

and [F2]fe,i = exp{~j2T:lk/N) (4) 
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where n, k — 0,1, . . . , N — I and 1 = 0,1,..., Lm — 1- Taking N samples of ry{n) in (|2]l, 

=D(e,77)Fi(7,)G((?)Xh„ +w,„, (5) 

where D(£,?7) = exp(j27re,,/A^), . . . , exp(j27re^(iV - l)/iV)], (6) 

G{0) =dmg[l,exp(-j27r6l/Af),...,exp(-j27r(7V-l)6l/7V)], (7) 

[huJjvAfTxi = [hf'^^, hf'^^, . . . , h^^^^]^, and [XjATxjvAfr = [^i; ^2, • ■ • , Xjv^] with 

X„ = diag[xu{0),Xu{l),...,Xu{N - 1)], and h„,„ = [ft.„^.u(0), /i„,„(l), . . . , ft.„,„(iV - 1)]^. We have CIR, 
denoted by [h„_„]L^xi, and CFR, denoted by [h„^]Arxi, related as h„,„ = F2h„ „. 

Thus, = D(e,7y)Fi(?7)G(0)X(I^v^ ® F2)h, +w, (8) 

where [h^jAr^i^xi = [hi'^^, hf^„, . . . , h^^ ^Y' . Stacking the outputs of Nr receive antennas, denoted by [r]Ar7v„xi = 
[r^, vl,..., r^^]^, and simplifying, 

r = A(e,77,6')h + w (9) 



where A(e, 77, 6) = I^^ ® (D(e, r;)Fi (77)G(0)X(I^^ ® F2)) (10) 



and [h\NnNTL„,xi = [hf , h^^, • • ■ , ^^^V ■ 

Neglecting the effect of SFO in the system model in (|9]l, i.e. letting ?; = 0, will result in the system model as 
given in and Q, where the effects of CFO, STE, and channel are considered. Similarly, if the STE is neglected 
in Q, i.e. making 6' = 0, we obtain the system model as in pl, where the effects of CFO, SFO, and channel are 
considered. Also, if the effect of CFO is not considered in (|9]l, i.e. making e = 0, will result in the system model as 
given in |8 |, for a SISO-OFDM system. Thus, the system model in (|9]l is general and considers the synchronization 
impairments and channel together 

III. CRLB Analysis 

A closed form expression of the CRLB for the estimation of CFO and channel is derived in IIT6II and 111 for 

the cases of single carrier and multicarrier communication systems, respectively. The related results can also be 

found in lfT4ll and ifTTll. Also, a joint estimation algorithm is described in lITSl . with CRLB derivation for the joint 

estimation of CFO, SFO, and channel for a SISO-OFDM system. But, calculation errors of equations (28) and 

(29) in IHl were reported and re-derivation of CRLB for the joint estimation of CFO and SFO for a SISO-OFDM 

system is done in ID, without considering the channel as a parameter to be estimated. In this section, we obtain the 

CRLB for the estimation of CFO (e) and SFO (rf) for a MIMO-OFDM system, considering the effect of STE (9), 

and analyzing the cases where channel is considered as a parameter to be estimated. The parameter vector of interest 

is represented as a = [e, rj, h]J, hj]^, where and h/ represent real and imaginary parts of h, respectively, and 

being a discrete parameter is omitted from the parameter vector of interest. Let Xi = X(IjVy ® F2), and /x be 

the mean of the received signal vector r in (|9]l. Then, the Fisher Information Matrix (F) is given by lfT6l . lfT9l . 

2 ^ \dii" da 1 
af„ oa. oa^ 
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where /x = (Iat^ (DFiGXi))h. 



(12) 



A. CRLB without channel being a parameter to be estimated 

Let Vuioc denote T without considering the channel as a parameter to be estimated. From ( fTTT i. 



Te,e Te.»7 

Let [Ciljvxiv = rfzag[0, 1, 2, . . . , (iV - 1)]. 



Then from ^ and (IT2] | we have 



(l + 77)DCi 



and 



9/1 

97 



j27r(l + r/) 

N 



DCiFiGXi h 



(13) 
(14) 

(15) 
(16) 



Substituting ( fT6] l in 7^ 
derivatives ll20l we get. 



9/x 

9e 9e 



and simplifying using the properties of Kronecker product and matrix 



(17) 



Let [C2]n-x.n = {diag(l]s) ® [diag{Ci)]^) o ([diag(lAr)]^ ® diag{Ci)) 
Then from Q, and ( fT2] i. 



9/1 

977 



9D 727re_^ , 9Fi j27r,^ ^ , 

where — = ^DCi, and —1 = :i_ C2 o Fi . 
orj N or] ly 



(18) 



Using ( fTSI ) and ( fTsT i, we obtain the closed form expressions for 7^ = — - — — — , 7^ ,: — 



drj drj 



de drj ' 



dr] de 



and 7,,,, 



as, 



7.,^ = 7,^e - ( I < 



^ Xf G^Ff Ci(eCiFi + (C2 o Fi))GXij j h (19) 



7-7 



^ = ® ( (^^^ Xf G^ (eFf Ci(C2 o Fi) + (C2 o Fi)^(C2 o Fi) + e^Ff C?Fi 

+ e(C2oFi)^CiFi)GXi)^h 



(20) 



Using ( flTl i. ( fT9l l. and (|20] |. the CRLBs for the estimation of e and r/, without channel being considered as a 
parameter to be estimated, denoted by CRLB(e^oc) and CRLB(77^oc), respectively, are obtained as. 



CRLB(e^oc) 
and CRLB (77 

woe J 



(21) 
(22) 
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B. CRLB with channel being a parameter to be estimated 

Let Tuic denote F considering the channel as a parameter to be estimated. From (fTTT) and ( fT3] l. 



From (fT2b we have, 



-5R 



lh,h 



(23) 



where Ih.h 



9/1 



9h, 



I®(DFiGXi), 



oh/ ohfl 



(24) 

(25) 
(26) 



Substituting (dSll, (|25ll, and ^ in Fh«,h« = -£—7^^ ^h,.h, = 1^1^' Th^.h^ = -£—7^^ ^h,.hR 



get. 



9h7 dhj' 



de dhl' 



Te.hi = tttt' = ^7—^' and simphfymg we 



dh.1 de 



IhH.hR = I ® (Xj G Fi FiGXi) Fh^^h^ = -jliiH.h^ = jIhj.hR, 



=h« I 



Xf G^Ff CiFiGXi ) ) = F^^^^ = = 



(27) 
(28) 



Similarly, using (fTSl l. ( l25T l. and (|26] | we obtain the closed form expressions for IJ,.h„ = 
I^,hi = — ^T— , and Fh^,^ — — — as. 



dfi 

drj dh.R ' 



dh.Fi d-q' 



drj dh.1 



dhj drj 



Ty^R.r, = ( I ® ( ( ^) Xf G^Ff (eCiFi + (C, o F^)) GX^) ) h = %1 = 1^,, = ^'^-^ 



(29) 



Using (fTsT l. dZTb . ( |28]) . and ( |29] l. the CRLBs for the estimation of e and 77, with channel being considered a 
parameter to be estimated, denoted by CRLB(e^c) and CRLB(?7^c), respectively, are obtained as, 

CRLB(e,„) = [F-i]i,i, (30) 



and CRLB(77^e) = ^'12,2. 



(31) 



Also, using MATLAB notation, we have CRLB for the estimation of h as. 



CRLB(h) - [F^c]3:(2L„+2),3:(2L„+2)- 



(32) 
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C. Significance of Joint Estimation 

The non-zero off-diagonal elements in (|23] | show the coupling between the impairments, CFO, SFO, STE, and 
channel, as observed in ([Til, which emphasizes the significance of a joint estimator. A study of coupling between 
CRLBs for the estimation of CFO and channel is explained in 1161 . without considering the effect of SFO and 
STE. In this section, we present coupling effects due to STE and channel, in the CRLBs for the estimation of e 
and 77 as shown in Fig[T] The CRLB equations in (EB, (ES, (O, and (EB are evaluated for a 2 x 2 MIMO-OFDM 
system with N — 128, and Lm = 6, having impairment values, e = 0.21, r/ ~ 120 ppm, and 6' = or —20, 
as shown in the figure. The derived CRLB expressions depend on the specific channel realization. Therefore, the 
CRLB expressions are numerically averaged over lO'^ channel realizations, using Spatial Channel Model (SCM) 
specified by the 3'''^ Generation Partnership Project (3GPP) ||2T1 . For the above given simulation setup, the change 



Coupling effect in a 2x2 MIIVIO-OFDIVI system 




10 12 14 16 18 20 



SNR{dB) 

Fig. 1. CRLB for the estimation of e and as a function of SNR(dB) sliowing the coupling effect. 

in CRLBs plotted in Fig[T] can be approximately represented as, 

SNR at CRLB(e^c)|e=o ~ SNR at CRLB(e^oc)|e=o + 6.25 dB 
SNR at CRLB(7y,„c)|0=o « SNR at CRLB (77^,0^)10=0 + 4.5 dB 
SNR at CRLB(e^oc)|e=-20 « SNR at CRLB(e^oc)|e=o + L25 dB 
SNR at CRLB(7y^oc)|e=-20 « SNR at CRLB (77^00)! 0=0 + 1 dB 

The above expressions show the effect of coupling of channel and STE on the estimation of CFO and SFO, 
respectively. Due to the coupling effect, the estimation of parameters without considering all impairments together 
results in performance degradation. 
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IV. Maximum Likelihood Estimation 

The ML cost function [3], |fT9l of the parameters e, 9, i], and h, obtained from (|9]l is, 

1 f - II r - Ah |P 1 
arg max P(r|e, 9, rj, h) = arg max - — exp <^ 5 ^ } . 

Computing the log-likelihood function and simplifying, we obtain an equivalent cost function, 

arg min J(e, 77, 0, h|r) = arg min (r — Ah)^(r — Ah). (33) 

The multi-dimensional minimization in ( |33T l gives the estimate of the parameters, e, 9, ij, and h, which is practically 
not feasible. With perfect channel knowledge at receiver, the optimization problem in (|33] l reduces to a three- 
dimensional minimization problem as argmin J(e, 77, 0|r, h), which also turns out to be too complex for practical 
purposes. In this section, we propose an ML algorithm in which the multi-dimensional optimization problem in 
(|33l l is reduced to a two-dimensional and a one-dimensional grid search, by rearranging the system model. Also, 
we propose a low complexity MML algorithm, which involves multiple one-dimensional searches, by making 
an approximation. Further, SML algorithm using different existing algorithms, which estimate fewer number of 
parameters, is also proposed. 



A. Proposed ML Algorithm 

The system model in ^ can be rewritten as 

r = (Ijv^ ® (DFiGX(Ijv^ F2))) h + w (34) 

= (Ijv, <E> (DFiGX(Ij^^ F2e„„J)) H + w, (35) 

where [Fae^Jfej = exp(-j27r/fc/7V), (36) 

and h = [hi^,l4^,...,H^J^, (37) 

with = [Hi%, U,^^, HJ^ J^, and h„ = [hl^ Oj^^^^^f for fc = 0, 1, . . . , iV - 1 and Z = 0, 1, . . . , (L + 
Omax — 1)- The system model in ( |35] | can also be represented as, 

r = (DFiX(IjY^ F2e„„J)) hg + w (38) 

= Aihe+w, (39) 

where Ai = Ijv„ (g, (DFiX(I^^ ® Fze^^J), (40) 

and he = [h?;,hf^,...,h?;^J^, (41) 

with h„, = [hf^_^)^,hf2^^)^,...,hf^^^^)J^, and h(„,„), = [Oj^^ hj,„ Ofg_^_g)^J^, 

for u = 0, 1, ... , Nt — 1, and v = 0,1,..., Nn — 1. From (|39] l, the least squares estimate of hg is given by, 

he ~ A\r. Let be the projection matrix of Ai. Therefore, the estimate of e and rj can be obtained from (3% 
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as. 



[tML,ffML\ = argmin(r - AiA|r)^(r - AiAjr) (42) 
= argmax || PA,r (43) 
= argmax Ji(e, r/|r) 

where Ji{e,ij\v) =|| PA,r \\l . (44) 

From (|9]), the least squares estimate of h is given by, h ~ A^r. Let be the projection matrix of A. Therefore, 
using the estimates of e and 77 from ( |43] ). the estimate of 9 can be obtained as, 

[Oml] = arg min(r - AA^r)^ (r - AA+r) (45) 



= argmax || I\r II2 (46) 
= argmax J2(6'|r, eML^VMh) 



where J2{B\r, eml, Vml) ^\~PAr \\l ■ (47) 
Finally, using the estimates of e, 77, and 9, we get the estimate of h as, 

[hAfi] = Atr (48) 

where A = Ijv« (D(eAfL, ??AfL)Fi(r)ML)G(^Mi)X (1^^ O) F2)) . (49) 
The series of steps involved are given in Algorithm 1. 

B. Proposed Modified ML (MML) Algorithm 

The estimation of e and rj, using ( l43T l. by a two dimensional grid search is not desirable for practical applications. 
Hence, a reduced complexity MML algorithm is obtained in this section. From (|40] l, we have 

Ai = Ijv^ ® (DFiX(I^^ F29„„J) 

= Ijv«®DA2, (50) 
where A2 - FiX(I^^ ^ F2e,„„J. (51) 

Therefore, from ( |38] l, (|40l i, and ( |50] l, the least squares estimate of hg is given by, 

he = A\r 

= (Af Ai)-iAf r 

= (Ijv, ® ((Af A2)-iAf D^)) r. (52) 
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Algorithm 1 ML Algorithm 



Inputs: N, Nt, Nr, r, X, L„ 



1: Evaluate F2, F2e„„^; 

2- for i ^rain • ^grid • ^vriax ^O 

3- for j 'Uniin ■ Vgrid • Vmax ^O 

4: Evaluate D(i, j), Fi(j); 

5: Ai(i, j) = Iw„ ® (D(i, j)Fi(j)X(IjVr »F2e_J) ; 

6: Ji(*,j|r) =11 PA,r Hi; 

7: end for 

8: end for 

9: [eA/Lj^L] = argmax Ji(j,j|r); 

10: Evaluate T>{€ml,t1ml), ^liVML); 

11: for fc = 0min '■ (^grid '■ (^max dO 

12: Evaluate G(fc), A{lML,f}ML,k); 

13: J2(fc|r, ^ML.flML) =|| 

14: end for 

15: [^a/l] 

16: Evaluate G{Oml); 

17: A = Ia,„ ® {p{€ML,TlML)^l{f]ML)G{OML)M^NT ® ^2 

18: [h.ML] = Atr; 

Output: [0ML , CML , flML , h^/i] 



arg max J2 (fc | r , eml , f]ML ) ; 



★ using (IHi, dlHl 



★ using (l6]l, (OJ 
* using (|40] | 
★ using (l44l i 



★ using (|6]l, (l3]l 

* using ([Toll 
★ using (|47] | 



★ using d?) 
★ using ( I49I 1 

★ using (l48T l 



Using (l50b and (|52] ). the cost function of the two dimensional grid search in ( |42] | can be written as, 

(r- AiAlr)^(r- AiAjr) 
= (r - ® (DA2(Af A2)-iAf D«)) r)"" (r - (l^, ® (DA2(Af As)"! A2«D«)) r) 
= r^(lAr^ (8) D) (Iwath - (Iat^ «) PaJ) (IwtVh - (Iwh ® PaJ) (Iw^ ® D^)r 
= d^R^ (Iatjvh - {Inu ® PaJ) - (Ijv„ PaJ) Rd* 

= d^CC^d* (53) 
where R = diag(r), d = diag{lN^ ® D), D^D = DD^ = Iw, Pa^ = A2(Af A2)"^Af , and 

C = R« {iNNn ~ (Ia^« ® Pa J ) . (54) 
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Algorithm 2 MML Algorithm 



Inputs: N, Nt, Nr, r, X, [0min,O^ax,Ogrrd], 

1: Evaluate F2, F2e„„^, Ci; 
2: R = diag{r); 

3- for j Vmin • Vgrid ■ Vriiax ^O 

4: Evaluate Fi(j), A2(j); 

5: C = (IjvAT^ - {Inj, ® Pa J); 



2tt(1 + j) 

6: ci(j) = -^^ dzagCIjvH ® Ci); 

_ cf^S(CCg)l 

8: Evaluate J3(j|r); 

9: end for 

10: fefMi] = argmin J3(j|r); 

j 

11: [^MMl] — ^{VMMl)', 

12: Evaluate Ti{eMML,TlMML), ^l{f]MML)\ 

13: Remaining steps are same as steps from 11 to 18 in Algorithm 1. 

Output: [9 MML , ^MML , fjMML , ^MMl] 



using dUl, (O, (Ell 



★ using (|3]l, (|5B 

* using (O 

* using dSST l 

* using dSTl l 
★ using 



★ using 
★ evaluation similar to step 7 
* using (l6]l, © 



For small values of e, d « 1 + jeci, where 

ci - ^^^^y^rfza.9(lAr^ ® Ci) (55) 
with Ci as given in ( fT4l i. Therefore, ( |53] ) can be written as, 

d^CC^d* « [1 + jeci]^CC^[l - jeci] 

= l^CC^l + jecf CC^l - jel^CC^Ci + e^cf CC^Ci 

= l^CC^l + e^cf 3fi(CC^)ci - 2ecf 3(CC^)1 (56) 



Differentiating ( |56] | with respect to e and equating to gives the estimate of e in terms of -q as, 

' ~ 5R(CC«)ci 

Substituting ( [57] l into ( [56] l, we have 



(57) 



J3(?7|r) = l^CC^l + e^cf 5R(CC^)ci - 2ecf 3(CC^)1. (58) 
Using ( fSSl ). we get the estimate of 77 by MML algorithm as 

[jyMAfL] = argmin J3(77|r). (59) 

»? 

Substituting ( [59] ) into (|57] | gives the estimate of e. Using the estimates of e and 77, the estimates of 8 and h can be 
obtained using (l47T l and ( |48] |. as done in Algorithm 1. The series of steps involved are given in Algorithm 2. 
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Remarks: For the range of values of |e| less than 0.10, the equation ( |56] l holds reasonably good with average 
error of the approximation less than lO^"'. 

C. Proposed Stage-wise ML (SML) Algorithm 

We also propose SML algorithm in which the joint estimation of CFO and STE is done in the first stage using 
the algorithm in [6i, ignoring the effect of SFO. Using the system model in (|9]l and ignoring rj, we have 

A(e, 0, 0) = In^ (g, (D(e, O)Fi(O)G(0)X(I^^ <g> F2)) (60) 
and [esML, Osml] = argmin(r - AAtr)^(r - AA^^r) 

€.6 

= argmax II PaF II2 (61) 
= argmax J4(e,6'|r) 

where J4(e,T?|r) Paf II2 . (62) 

With the estimates of CFO and STE from ( |6TT i in the first stage, the received signal is used to estimate SFO and 
channel using the ML algorithm in [8|, as extended to a MIMO-OFDM system, in the second stage. 

[flSML] = argmin(r - AAtr)^(r - AA^r) 

n 

= argmax || Paf \\l (63) 

= argmax J5(r/|r, isML,QsML) 

where Jsivl^, esML, Osml) ^\\ \\l . (64) 

Finally, using the estimates of e, rj, and 9, we get the estimate of h as in ( l48T l. The series of steps involved are 
given in Algorithm 3. 

D. Comparison of Computational Complexity 

The computational complexity of ML algorithm depends mainly on computation of the two-dimensional grid 
search in (l43T l. for obtaining the estimates of e and rj, and the one-dimensional grid search in ( |46] |. for obtaining 
the estimate of 6. Let g^, g,^, and gg denote the number of grid points used in search for e, 77, and 6. Therefore, the 
computational complexity of evaluating ( |43] l is approximately equal to g^gj^O{N^), whereas that of evaluating ( |46] l 
is approximately equal to ggO{N^). Thus, the total computational complexity of ML algorithm is approximately 
given by {g,g^ + ge)0{N^). 

The computational complexity of MML algorithm depends mainly on the computation of two one-dimensional 
grid searches, given in ( |59] l and (l46T l. The computational complexity of evaluating ( |59] l is approximately equal 
to gfjO{N'^), whereas that of evaluating (|46] | is approximately equal to ggO{N^). Thus, the total computational 
complexity of MML algorithm is approximately given by (g,, + gg)0{N^). 
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Algorithm 3 SML Algorithm 



Inputs. A^, A^Ti A^_R? r, X, L^a^ [(^min7^max^^grid\i {^min: ^max: ^grid\i {'Hmin : "^max ^ Vgrid] 



1: Evaluate Fa, Fi(0); 

2: for i — ^rnin • ^grid • ^riiax dO 



3: for k 



^ arid ■ ^rnax dO 



4: Evaluate D(i,0), G(fc); 

5: A(z, 0, k) = ® (D(z, 0)Fi(0)G(fc)X(I,v^ ® F2)) 

6: J4(i,fc|r) =11 PAr ||2; 

7: end for 

8: end for 

9: [esAfL,^SML] = arg max J4(i, fc|r); 

i,k 

10: Evaluate G(^sAfL); 

11: for j = rjrnin ■ Vgrid ■ Vmax dO 

12: Evaluate D(esML, j), Fi(j), A.{esML-,i,OsML); 
13: J5U\r,esML,0sML) =\\ 
14: end for 

15: [vsml] = argmax J5(j|r, esMLi^'sA/i); 

16: A{isML,VSML,dsML)', 

17: [hsML] = A^r; 

Output: [6sML , ^SML ,flSML, hgMi] 



7k- using (|4|i, (|3]l and rj is ignored 



★ using (|6]l, (|7]i and 77 is ignored 
★ using (|60] | 
★ using 



★ using d?) 

★ using dg, (Oil, (HOll 

★ using (|64] | 



★ using (fTOl l 
★ using 



The computational complexity of SML algorithm depends mainly on the computation of the two-dimensional 
grid search in dMT l. for the estimation of e and 9 in the first stage, and the one-dimensional grid search in ( |63]) . for 
the estimation of rj in the second stage. The total computational complexity of the first stage is approximately equal 
to gegeO{N^), whereas that of second stage is approximately equal to griO{N^). Thus, the total computational 
complexity of SML algorithm is approximately given by (g^ge + gri)0{N^). 

Remarks: Based on the system setup in Section V, the typical values of g^, g^j, and gg are 81, 101, and 
43, respectively. Thus, MML algorithm is around 57 times and 25 times faster than ML and SML algorithms, 
respectively. 

V. Simulation Results and Discussions 

A. System Setup 

The simulated 2x2 MIMO-OFDM system has — 128 subcarriers for each transmitter with 20 MHz signal 
bandwidth. The SCM specified by 3GPP ETI . is used to generate the fading channels in the Urban Macro Scenario 
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with L„i = 10. Also, The transmitted symbols belong to QPSK constellation with unit amplitude. Therefore, the 
variance of the complex additive Gaussian noise at the receiver, ct^, is varied for getting SNR(dB) in steps of 5 for 
the evaluation of Mean Square Error (MSE) and CRLB. Also, lO'^ trials are carried out in simulations for getting 
each point in the MSE and CRLB plots. 

We consider the training blocks having a CP of length 32. The condition (L^ + Omax) less than length of 
CP results in Omax 21 and 16*1 < 21. The range of normalized CFO used for grid search is |e| < 0.4 with a 
resolution of 10^^ and that of normalized SFO is jT^j < 5 x lO^'^ with a resolution of 10^^. The actual values of 
the impairments, e, r/, and 6 used in the simulations are 0.021, 0.000101 or 101 ppm, and 2, respectively. 

B. Performance Measures 

The ML estimates of the parameters are used for calculating the MSE values as, 

Ntrials 

E llp.-plli 

MSE(p) = . (65) 

trials 

where p represents the actual parameter, p^ represents the estimate of the parameter at i*'^ trial, and Ntrials represents 
the number of trials. The Probability of Timing Failure [Ptf) is used as an indicator to illustrate the robustness of 
symbol timing estimator Q, which is expressed as. 



Ptf{p)=Pr 



-9 :^ 



(66) 



where p is the absolute difference between the estimated value of 9 and the actual value of i 



C. Performance Assessment 

The MSE values of the estimated parameters are calculated using i65[ and are plotted in log-scale against 
SNR(dB) for ML, MML, and SML algorithms in Fig.2-Fig.5. The CRLBs of the parameters are also plotted in the 
corresponding figures. It is found from Fig|2]and Figl3]that the MSE plots of the ML algorithm for the estimation 
of 1] and e closely follows CRLB(?7^c) and CRLB(e^c), but with a performance degradation of around 1.25 dB 
SNR and 0.75 dB SNR, respectively. Also, it is found that there is only small performance difference between ML 
and MML algorithms. Performance loss of more than 3dB, and error floor at 25 dB occurs for SML algorithm 
when compared with the performance of ML or MML algorithms in Fig|2] and Fig[3] respectively. 

The probability of timing failure for the estimation of 6 is calculated for p — 1, as defined in (|66] |. and is 
plotted in Fig|4]for ML, MML, and SML algorithms, respectively. As in the cases of e and ?/, there is only small 
performance degradation due to the approximation in the MML algorithm, and performance loss of more than 3dB 
occurs for SML algorithm when compared with ML algorithm, as observed from the plots in FiglH The MSE plots 
of ML estimates of h using ML algorithm with STE=0, denoted as ML(6' = 0), together with ML, MML, and SML 
algorithms are shown in Fig|5j Tr(CRLB(h)) evaluated using ( |32] | is also shown in the same figure. It is found from 
Fig|5]that MSE plot of ML(6l = 0) follows Tr(CRLB(h)) with a performance loss of 4 dB. Comparing the MSE 
plots of ML(6' = 0) and ML in the figure, it is found that the performance degradation of ML algorithm below 20 
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Synchronization lmpairments-(e=0.021 , ti=101 ppm, 9=2) 
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SNR(dB) 



Fig. 2. CRLB and MSB for the estimation of SFO as a function of SNR(dB) using ML, MML, and SML algorithms. 



Synchronization lmpairments-(e=0.021 , ri=101 ppm, 6=2) 
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Fig. 3. CRLB and MSB for the estimation of CFO as a function of SNR(dB) using ML, MML, and SML algorithms. 



dB SNR is due to the coupling effect of STE on the channel. As in the cases of e, ?; and 9, there is only small 
performance degradation due to the approximation in the MML algorithm. Also, the performance degradation and 
error floor of the MSE plot of SML algorithm is due to the stage-wise estimation of the parameters. 

Thus, it is found from the figures that there is only small performance degradation due to the approximation in 
the MML algorithm when compared with ML algorithm. Also, performance degradation occurs, if joint estimation 
of all impairments together is not done, as observed from the performance plots of SML algorithm in the figures. 
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Fig. 4. Probability of Timing Failure as a function of SNR(dB) using ML, MML, and SML algorithms. 



Synchronization lmpairments-(e=0.021 , ri=101 ppm, 9=2) 




5 10 15 20 25 30 35 



SNR(dB) 



Fig. 5. CRLB and MSE for the estimation of channel as a function of SNR(dB) using ML, MML, and SML algorithms. 



VI. Conclusion 

In this paper, ML algorithms, for the joint estimation of synchronization impairments and channel in a MIMO- 
OFDM system, are proposed. A system model, which shows the effects of the synchronization impairments and 
channel, is formulated. CRLBs for the estimation of continuous parameters are derived, which show the coupling 
effect among different parameters and the significance of a joint estimator. An ML estimation algorithm is proposed, 
where the multi-dimensional optimization problem for estimating the parameters is reduced to a two-dimensional 
and a one-dimensional grid search. To reduce the complexity further, MML algorithm which involves only multiple 
one-dimensional searches is also proposed. Further, SML algorithm using existing algorithms, which estimate fewer 
number of parameters, is also proposed. The performances of the estimation methods are studied through numerical 
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simulations and it is observed that, there is only small performance degradation due to the approximation in 
MML algorithm. Also, the proposed ML and MML algorithms exhibit better performance than the proposed SML 
algorithm, which uses existing algorithms. 
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